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The measurement of the top forward-backward asymmetry in tt production measured at the 
Tevatron shows deviation from the standard model prediction. A u — > t transition via a flavor- 
changing Z' can explain the data. We show that left-handed t^u^Z' couplings can be constrained 
from Bd,s mixing while the constrains on the right-handed couplings tRUaZ' vanish in the limit of 
m u — > 0. We then consider the most general form of the tuZ' interaction which includes vector-axial 
vector as well as tensor type couplings and study how these couplings affect the top forward-backward 



m 

asymmetry. 
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\ The top quark with its high mass may play a crucial role in electroweak symmetry breaking. Hence the top 

sector may be sensitive to new physics (NP) effects that could be revealed through careful measurements of 
top quark properties. The top quark pair production in proton- antiproton collisions at the Tevatron collider 
with a center-of-mass energy of ^/s — 1.96 TeV is dominated by the partonic process qq — ?> tt. Recently 
the CDF experiment has reported a measurement of forward-backward asymmetry in tt production which 
appears to deviate from the standard model (SM) predictions. The CDF collaboration measured the forward- 
backward asymmetry (Afb) in top quark pair production in the tt rest frame to be A l p B — 0.475 ± 0.774 
for M t i > 450 GeV l|, which is 3.4 a deviations from the next-to leading order (NLO) SM prediction 
A l ^ B = 0.088 ± 0.013 JQ- The D0 collaboration also observed a larger than predicted asymmetry 0]. 
The current measurement of the top quark pair production cross section from 4.6 fb" 1 of data at CDF is 

u t i = (7.50±0.48)p6, (1) 



for m t = 172.5 GeV fl, in good agreement with their SM predictions by Langenfeld et al. a tt = 7.4Sit\[ 
pb Cacciari et al. a t i = 7.26lg'gg pb 0, Kidonakis a t i = 7.29lggg pb 10 1, and recent Ahrens et al 

10.31 



+0.66 
80 



significantly low value a t i = 6.30±0.19_q 2 3 pb [ll(- Hence new physics models that aim to explain the Afb 



measurement must not change the production cross section appreciably. Many NP models that affect Ap B , 
either via s-channel 0-0 

or i-channel exchange of new particles [3ll457j have been proposed to explain 
the forward-backward anomaly. Here we will focus on the model with a Z' boson that has a flavor-changing 
tuZ' coupling. This coupling can contribute to tt production at the Tevatron via the t-channel exchange of 
the Z' boson (see Fig. Ufa)). The Afb measurement can be explained with a light Z' with a mass around 
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FIG. 1: Left panel(a): Tree-level ti production diagram involving the Z' exchange. Right panel(b): Tree-level 
diagram with tq'Z' coupling (q' — u, c, t) which generates an effective bqZ' ( q =d,s) coupling through a vertex 
correction involving the W exchange. 



Flavor-changing neutral current (FCNC) effects in the SM are tiny and satisfy the condition of natural 
flavor conservation proposed by Glashow, Weinberg and Paschos [59(. The condition of natural flavor 
conservation can be avoided if quarks of the same charge couple to more than one Higgs or their couplings to 
a new vector boson (e.g. a Z' boson) are different for different generation. To date there is no experimental 
evidence of FCNC effects beyond those expected from the SM. There are some anomalies in the B system 
which might require new physics to resolve, but the NP-generated FCNC effects that are needed in the B 
system are much smaller than the one needed to resolve the top Afb [sB]. A tree-level dbZ' or a sbZ' 
coupling is strongly suppressed by Bd, s mixing. A tree-level tq' Z' coupling, where q' = u, c, t, will generate 
an effective bqZ'(q = d, s) coupling through a vertex correction involving the W exchange [61( (see Fig.[TJb)). 
The B q mixing constraints on these effective vertices would then lead to constraints on the tq'Z' coupling. 
The vertex corrections are divergent and can be regulated by a cut off A, which represents the scale of NP in 
an effective theory framework. In NP models where there are no bare bqZ' couplings, the vertex corrections 
with a chosen A can be used to constrain the tq'Z' coupling from B q mixing measurements. We will take 
the scale of new physics to be ~ TeV . In specific complete models A will represent the mass of some new 
particles. In models of NP where there are bare bqZ' couplings the vertex correction will renormalize the 
bare bqZ' vertices to produce the renormalized vertices U q b- These renormalized vertices can then be fitted 
to B q mixing data. Assuming the vertex corrections to be less than or at most the same size of the bare 
couplings one, we can obtain bounds on the tqZ' couplings by requiring the generated bqZ' coupling to be 
< U q b- It is possible to have models where large bare bqZ' couplings cancel with large vertex corrections 
to produce small renormalized bqZ' vertices consistent with experiments. We will not consider these finely 
tuned model. 

When the vertex corrections are computed, one finds that right-handed tuZ' couplings do not contribute 
to B q mixing in the limit of setting the up quark mass to zero. We note that ttZ' couplings do not have such 
suppression and will contribute to B q mixing via the vertex corrections. Even though the ttZ' coupling does 
not contribute to the top Afb, in specific models of NP this coupling may be related to the FCNC coupling 
tuZ' 62j . It turns out the B q mixing constraints on ttZ' are weak because of the small Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements V ts (d) and not because of right-handed couplings. The tqZ'(q — u,c,t) 
couplings via box diagrams can produce an effective d(s)buu operator that can contribute to decays like 
B — s- K(K*)ir(rj,ri' p) or B — > n(p)ir(p), etc. decays. The effects of these new operators can be observed in 
CP-violating and/or triple product measurements [63[. However, these effective operators only modify the 
SM Wilson's coefficients in the SM effective Hamiltonian and so the CP-violating predictions and/or triple 



3 



product measurements should be similar to the SM for a reasonable choice of tqZ'(q = u, c, t) couplings. 

We will next consider the most general tuZ' couplings including both vector, axial vector and tensor 
couplings ( ~ -7^7-) an d study the effect of these couplings on the top Aps- The interesting feature about 
these tensor couplings are that we can avoid the B q mixing constraints due to the suppressions of these 
operators at low energies j&i ]. The momentum dependence of these operators imply that at the b quark 
scale these operators will be suppressed by ~ m&/mt and consequently the B q mixing constraints will be 
weak for these operators. 

The paper is organized in the following manner. In the next section we discuss the B q (q = d, s) constraints 
on the tuZ' operators. In the following section we introduce the general tuZ'coupling including tensor terms 
and study the effects in the top Aps- This is followed by the section on the t — > uZ' branching ratio 
calculations. In the final section we present our conclusions. 



I. CONSTRAINTS ON tq'(= u,t)Z' COUPLINGS FROM B qi=d>s) MIXING 

In general, new physics contributions to the mass difference between neutral B q meson mass eigenstates 
(AM q ) can be constrained by the AM q experimental results. In the SM, B q -B q mixing occurs at the one- 
loop level by the flavor-changing weak interaction box diagrams. The mixing amplitude M* 2 is related to 
the mass difference AM q via AM q = 2|M' 2 |. The recent theoretical estimations for the mass differences of 
B®-B® and B®-B® mixing [(35| at la confidence level are 

(AM S ) SM = le.St^ps" 1 , (AM d ) SM = 0.555±° ;^ps- 1 . (2) 

The latest measurements of mass difference by CDF [66| and D0 [67[ for B s mixing are 

AM Bs = (17.77 ±0.10(stat.)±0.07(syst.))ps~ 1 

AM Bs = (18.53 ±0.93(stat.)±0.30(syst.))ps^ 1 . (3) 

The Heavy Flavor Averaging Group value for the mass difference of B®-B® mixing is AAfs d (exp) = (0.507 ± 
0.004) ps _1 [sHl. The experimental results for the mass differences of both B®-B® and B^-B^ mixing are 
consistent with their SM expectations. Hence, the mass difference results can provide strong constraints on 
NP contributions. 

In this section we will consider the -B^.s mixing constraints on the tq'(= u,t)Z' couplings. 



1. tuZ' left-handed coupling 



The most general Lagrangian for flavor-changing tuZ' transition is 



C 



tuZ' 



m t 



(4) 



where q = pt — p u - In general, the couplings a, b, c and d are complex and can be momentum-dependent 
(form factors). In this work we will take the couplings to be constants with no momentum dependence. 
Consider the tuZ' vertex with a = —b = g^ u , and c=d=0 in Eq. ((H). This generates effective bqZ'{q = d, s) 
coupling at one-loop level due to W exchange. We obtain the bqZ' coupling in the Pauli-Villars regularization 
as 



Cz> =U qb qr(l~l5)bZ^ (5) 
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where 



U, 



qb 



9h 



V2 



x t Log[ 



(6) 



where A ~ TeV is a cutoff scale, and x t = mf/M^- The function U q b includes only the contribution from 
the W boson, and the contribution of the associated Goldstone boson in the SM is the order of m u /Mw- 
Note that for Bd mixing the coupling g^ u is associated with the CKM factor V* d Vtb ~ 1, and thus one can 
expect a strong constraint on g^ u from the mass difference AMd- 

A tree-level exchange of the Z' generates the AB = 2 effective Lagrangian responsible for the neutral B q 
meson mixing 



(7) 



where (qb)v-A = 37^(1 — l5)b, and the QCD correction factor r\z> = [a s (Mz')/a s (m b )] e ^ 23 . The Z' 
contribution to the B q mixing amplitude can be obtained by using the vacuum insertion method as 



M 2 f = \^z. mBq f B B q . 



(8) 



In the presence of new physics, the mixing amplitude M\ 2 can be parameterized by complex parameters 



M? 2 = [MU SM A q . (9) 

In our case, A q = \A q \e 1 ^ = 1 + \M? ] Z ' /\M? ] SM . A global analysis on the parameters |A g | and <t> q for 
Bd — Bd and B s — B s mixing are carried out in [65j . The best fit results for Ad and A s in this analysis at 1 
a confidence level (scenario I) are 

|A d | = 0.747t° ;^ 5 2 , ^ = -12.9± 3 2 ;f , (10) 



and 



A S |= 0.887±°;S& ^ = -51.6tft a * or -130.01^. (11) 



The A d constraint in Eq. {TDJ on the coupling g^ u at m t (m t ) = (165.017 ± 1.156 ± 0.11) GeV [65 1 



pSM _ 27.2t\\ [65J, and M z > = 150 GeV is shown in Fig. [2J The numerical values of all other theoretical 



inputs can be found in 65j. They are varied within 1 a errors in the fit. The cutoff scale A is varied between 
300 GeV to 2 TeV . The green scatter points in Fig. [2] satisfy only |Ad| in Eq. (fT7J|) . while blue points satisfy 
both | Ad | and 4>f in Eq. flTUl) . The results indicates that Bd mixing can strongly constrain the tuZ' coupling 
g^ u even at A = 300 GeV. In particular we note that the maximum value for \g^ u \ is around 0.2 and is 
associated with a large phase. In fact there are no real g^ u that satisfy the Bd constraint. 

On the other hand, Fig. [3] suggests that the constraints from B s mixing on the tuZ' coupling g^ u are 
weaker (~ 0(1)) even at A = 2 TeV. This can be understood from the fact that the B s mixing contribution 
in this case is associated with the CKM factor V* s Vtb and is suppressed. The (green, blue, red) scatter 
points in Fig.[3]are constrained by ( |A S |, {|A S |, 4>f = -51.6±g 4 7 2 °}, {|A S |, </£ = -130.0±^°}) in Eq. (fTT|) . 
respectively. The large negative phase tfif' prefers large g^ u values. 
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For Ad constraint 




Arg[g^][Deg] A [GeV] 

FIG. 2: \gtu\ vs Arg[g^] [Deg] (left panel) and \gtu\ vs A[GeV] (right panel) for Bd mixing. Green scatter points are 
constrained by |Ad|. Blue scatter points are constrained by |A<j| and fa , 



For A s constraint 




FIG. 3: \gtu\ vs Arg[<7^,] [Deg] (left panel) and |p tu | vs A[GeV] (right panel) for B s mixing. Green scatter points 
are constrained by |A S |. Blue scatter points are constrained by |Ad| and </>^ = — 51.6lg 4 7 2 . Red scatter points are 
constrained by |A<j| and <j>% = — 130.0lj2 • 



2. tuZ' right-handed coupling 

We now consider the tuZ' vertex with right-handed couplings, a = b = gji , and c=d=0. The contribution 
of this vertex to Myi is suppressed by m^/m^. Hence, the right-handed coupling gn cannot be constrained 
by B q mixing. 

Finally as indicated in the earlier section, the left- and the right-handed couplings generate via the box 
diagram effective qbuu{ q — d, s) operators. These operator can be constrained by observables in nonleptonic 
B meson decays like B —> tttt/Ktt. These operators change the Wilson's coefficients of the SM effective 
Hamiltonian with the change being ~ 10~ 2 at the scale fi — Mw for Mz> = 150 GeV and gt u >9tt ~ O(g). 
Since the generated NP physics operator structures are similar to the SM there are no easy way to detect 
their presence. A detailed fit to all the nonleptonic data may provide constraints on the couplings Qp,', R , 
which we do not perform in this work. Some analysis along this line has been done for tdW 1 coupling in [70J • 



3. ttZ' coupling 

For completeness, next we consider B q mixing constraints on the ttZ 1 couplings. The Lagrangian for the 
ttZ' interaction is 

L ttZ ' = % t V(l-75)+^7 Al (l+75)]^;. (12) 



6 



Again, we evaluate the one- loop diagram (see Fig. [Tf b)) in the Pauli-Villars regularization and obtain the 
effective Lagrangian for bq(=d,s)Z' interaction as 



L, z , 



(13) 



where 



U' qb 



V2 



with 



(4 



1 f 1 , f l ~ V ,\L( T r XA 2 , 1 X 2 \ R (\ T . xA 2 , 1 

^)Jo dx J dy i gtt \ Log[ K^ t ] + 2 d7J + ^lo 9 [ T7 ^] + 



D, 



(14) 



(15) 



and D u = x + (1 — x)xt- The function f tt includes both the W boson and the associated Goldstone boson 
contributions. The ttZ' contribution to the B q mixing amplitude is 



[M- 



q ]Z' 



A[U '" b? n m f 2 B 



(16) 



Both Bd — Bd and B s — B s constraints in Eqs. (fT0|) and (fTTj) can allow large ~ 0(1) values for g tt ' . 



II. TOP QUARK FORWARD-BACKWARD ASYMMETRY 



In this section we calculate the top Aps keeping in mind the constraints derived on the coupling from 
the previous section. The most general Lagrangian for a flavor-changing tuZ' interaction is given in Eq. 
This interaction can contribute to uu — > tt scattering amplitude through the t-channel exchange of the Z 1 
boson (see Fig. Ufa)). The tree-level differential cross section for qq — > tt process in the ti center-of-mass 
frame including both the SM and for Z' contributions is 



da 
dcos 9 



327TS 



(A 



SM 



Asm-z 1 + Az 1 ) ■ 



(17) 



where s = (p q + p q ) 2 is the squared center-of-mass energy of the tt system, f3 t = y/l — 4mf/s, and the 
polar angle 9 is the relative angle between direction of motion of the outgoing top quark and the incoming 
q quark. The quantities Asm, Asm-z 1 , and Az' denote the leading order SM, the interference between the 
SM and Z\ and the pure Z' scattering amplitudes, respectively. These amplitudes can be obtained in terms 
of kinematic variables 9 and s as 



Asm 


_ 2 9t 
9 


AsM-Z' 


9 


A z , 


1 r 

~~ 4 Lf 



4mr 



Ml 



M 2 ,) 2 + M 2 ,r 2 z , 



(h + h) , 
{fa + h + h). 



(18) 



Where eg = (3 t cos9, and i = (p q — p t ) 2 = — 3/2(1 — f3 t cos9) 



The functions /jS (i = 1-5) can be found 



in the Appendix. Here we assume the couplings a, b, c and d to be real. Our results for tt production 
are obtained by the convolution of the analytic differential cross section of Eq. (IT71) with the CTEQ-5L 
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parton distribution functions [71] implemented in Mathematica. We expect the MSTW 2008 [72| parton 
distributions to give compatible results. 

The forward-backward asymmetry of the top quark in the ti cm. frame is defined as fl^ 



Att _ 

a fb — 



of — ob 



where 



a F = 



do 
o d cos i 



■dcos9 , <jb = 



da 



i dcost 



■dcosO. 



(19) 



(20) 



In our analysis, we choose some representative values for the couplings a, b, c, and d to generate large 
forward-backward asymmetry Ap B for high M t \ (> 450 GeV) without distorting the shape of the mass 
spectrum da t ij dM t i- We fix the renormalization and factorization scales at fiR = fiF — mt. We evaluate 
Ap B which includes the NLO SM and the Z' contributions at mt = 172.5 GeV. Also, we apply a QCD 
K- factor K = 1.3 to the tree-level cross section in order to match the SM prediction for a t i- We consider the 
Z' boson with mass Mz' = 150 GeV and width Tz> = for the numerical analysis. 



A. Pure vector-axial vector couplings: a = ^6 and c = d = 



This case has already been considered before [3l| , but only right-handed couplings were considered. Here 
we will consider both right- and left-handed couplings. We take the representative values of the couplings 
a = —b = \gt u \ = 0.257, and c = d = 0. This value for g^ M satisfies the |A<j| constraint but not the phase 
4>^ constraints from Bj, mixing (see Fig. [5]). For these values A l p B can be explained within one a error of 
its measurement for M t i > 450 GeV. In Fig. 01 we show the M t i distribution for the ti observables Ap B , 
and <r t t. The differential distribution, da t i/ dM t t, has been measured in eight different energy bins of M t j for 



m t = 175 GeV in Ref. 74]. Our distribution of da t i/dM t i is consistent with the measurements. Since the 



partonic scattering amplitudes in this case (see the Appendix) depends on b 2 and b 4 terms, our results hold 
for right-handed couplings also, i.e a = b = \gf u \ = 0.257, and c = d = 0. 



-b= \gt\= .257, c=d=0 



-b= \g^J= .257, c=d=0 



400 500 600 700 800 900 
M, r [GeV] 




400 600 800 1000 1200 1400 
M, F [GeV] 



FIG. 4: Left panel: M t i distribution of Ap B in the two energy ranges [350,450]GeV and [450,900]GeV of invariant 
mass Mu- Green band: the SM prediction. Blue band with la error bars: the unfolded CDF measurement [lj]. 
Red line: the SM with Z' exchange prediction for (a = — b — 0.257, c = d = 0). Right panel: M t j distribution of 
da t t/dM t t [in fb/GeV] for eight different energy bins of Mu- Green line: the NLO SM prediction. Blue band with 
la error bars: the unfolded CDF measurement Red line: the SM with Z' exchange prediction for above values 
of couplings at mt = 175 GeV. 
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B. General case: all couplings are present 

In this section we consider the most general tuZ 1 couplings. We showed earlier that the left-handed 
coupling are strongly constrained from Bd mixing and there are no real values of g\ u that satisfy the Bd 
mixing constraint. We now investigate the effect of the couplings c and d on the Afb predictions. 



C. Pure tensor couplings : a = b = and c = ±d 

We consider the case of pure tensor couplings. In this scenario we can avoid the B q mixing constraints as 
the effects of the tensor couplings are suppressed by ^ at the b mass scale. The SM and Z' interference 
contribution Asm-Z' in Eq. (Ti"%l) vanishes in this case. The functions f± and /s in pure Z' contribution Az 1 
are also zero, and f% is order of (cs/m t ) 2 . The mass spectrum for A l p B is shown in Fig. [5{a) for only c = ztd 
couplings (c = ±d — 0.5). The results indicate that Z' contribution cannot reproduce the Afb measurement 
within one a for M t i > 450 GeV even at a low Mz> = 100 GeV (yellow lines) value. 

a= b= 0, c=± d=0.S 



0.4 
0.2 



500 600 700 800 

M t ,-[GeV] 



FIG. 5: M t t distribution of A FB . Green band: The SM prediction. Blue band with la error bars : CDF measurement. 
Red and yellow lines: The SM with Z 1 exchange prediction at M Z i — 150 GeV, and M z i = 100 GeV, respectively, 
for a — b = and c = ±d = 0.5. 



D. All the couplings are same order 

Finally, we consider the case where all couplings are of the same order. We choose the representative 
values of the couplings a = — b = \g^ u \ = 0.239, and c = d = 0.148. Again this value for g^ u satisfies the 
| Ad | constraint but not the phase (j)^ constraints from Bd mixing (see Fig. [2]). In Fig. [BJ we show the M t i 
distribution for the tt observables A f p B , and cr t j . We note that A f p B can be explained within one a error of 
its measurement for M t i > 450 GeV. The distribution da t t/dM t i is also consistent with the measurements. 
Similar results are obtained with a = b = \g^ u \ = 0.245, and c = d = 0.148 as shown in Fig. [7J The 
conclusion is that the inclusion of the tensor couplings does not have a significant effect on the top Afb and 
can only slightly lower the values of the couplings a and b relative to their values in the pure case, with no 
tensor couplings, discussed earlier. The presence of the tensor couplings may have an important impact on 
the polarization measurement in ti production [7o1 |. 
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a=-b = 0.239, c=d=0.14 



400 500 600 700 800 900 
M, ,[GeV] 



a=-b = 0.239, c=d=0.148 



10 1 



- 10" 1 



10~ 2 



400 600 800 1000 1200 1400 
M, ,[GeV] 



FIG. 6: M t j distributions of A" s and da t t/dM t t [in fb/GeV]. Pink lines: the SM with Z' exchange prediction for 
(a — ~b = 0.239, c = d = 0.148). The same conventions as in Fig. 0]used for other lines. 



a=b = .246, c=d=0 .H 



400 500 600 700 800 900 
M, ,[GeV] 
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a=b = .246, c=d=0 .148 



400 600 800 1000 1200 1400 
M, ,[GeV] 



FIG. 7: M t t distributions of A" s and dcru/dMu [in fb/GeV]. Pink lines: the SM with Z 1 exchange prediction for 
(a — b = 0.245, and c = d = 0.148). The same conventions as in Fig. |4]used for other lines. 



III. 



t -> uZ' BRANCHING RATIO 



In this section we consider the decay width for t 
coupling is given as, 



r(f uZ') 



i 



lGirmt 



1 



uZ' . The decay width with the most general tuZ' 



(m 2 t +2m 2 z ,){a 2 +b 2 ) 



-6m 2 z ,{ac - bd) + m%,(—%- + 2)(c 2 + d 2 ) 



(21) 



Branching ratio is defined as 



BR 



tuZ 1 



r[t-> cZ'} 

T[m t ] 



For the top width we use T(m t ) ~ T(t — > bW) which is given by, 



T(t -> bW) = 



Gi 



8tt%/2 



\Vtb\ 2 m 3 t 



1 



m 



w 



1 



- 2- 



nv 



(22) 



(23) 



In Fig. [5] we show the variation of t — > uZ' branching ratio with Mz> for different couplings. For couplings 
a = ±b = 0.257, and c = d = (red dashed line), we get BR tu z' ~ 6% at m t = 172.5 GeV, for a = - b = 
0.239, c = -d = 0.148 (blue dashed line), BR tu z> is 6.9%, and for a = b = 0.246^ = d = 0.148 (pink dashed 
line), BRt u z' is 7.2%. These branching ratios may be observable at the LHC |58j . 
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FIG. 8: BR Z , vs M z >- Red dashed line is for a = ±6 = 0.257, and c = d =0. Blue dashed line is for a = -b = 0.239, 
c = -d = 0.148. Pink dashed line is for a = b = 0.246, c = d = 0.148. 



A large forward-backward asymmetry in tt production, about a 3.4<r away from the SM prediction, has 
been reported by the CDF collaboration. A Z' with flavor-changing tuZ' coupling can explain this anomaly. 
In this work we considered B^ tS constraints on the tq'Z' couplings ( q' = u,t). These constraints resulted 
from the bounds on the effective b(s,d)Z' vertices generated from vertex corrections involving the tuZ' 
couplings. We found that the right-handed couplings were generally not tightly constrained but the left- 
handed couplings were tightly bound from the B^ s mixing data. We then considered the most general tuZ' 
coupling including tensor terms and found that the tensor terms did not affect the top Apg in a significant 
manner. Finally we computed the branching ration for the t — > uZ' transition and found it to be in the 
percentage range. 

Acknowledgments: This work was financially supported by the US-Egypt Joint Board on Scientific and 
Technological Co-operation award (Project ID: 1855) administered by the US Department of Agriculture. 



For the scattering amplitudes calculation in tt center-of-mass frame, we choose the relevant coordinates 
of particle momenta as 



With this choice and assume all the couplings in Eq. Q to be real, we obtain the functions /, in the 
scattering amplitude in Eq. (|Tg)) as 



IV. CONCLUSION 



Appendix A: Functions in scattering amplitudes 



Pg,g = ~(1,0,0,±1), 

m = |(l,±Asin0,O,±/3tCos0). 



(Al) 



fi 





(A2) 
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h = T&i 



32(a 4 + 6 4 ) ( 3 + 2c e + c 2 g + 2ft 2 ) + —r (c 4 + d 4 ) (32(9 - 2c + c 2 g )m 



+32(-5 + 3c e + c 2 g + c 3 e )m 2 s + s 2 (5 - 2c 2 e + ft 2 - c e {3 + ft 2 )) 2 ) + 128a 3 c( - 2c e + c 2 , + ft 2 ) 
— ^2c 2 d 2 ( - 32(-5 + 3c e + c 2 + c 3 e )m 2 s + 32to 4 (-11 + 6c e + c 2 . - 4ft 2 ) 

-s 2 (5 - 2c 2 + ft 2 - c e {3 + ft 2 )) 2 ) + ^ac(8b 2 m 2 (-2 + 2c g + 3c 2 , - 3ft 2 ) 



5ci 



+c 2 (-l + c e )(8(-3 + c e )m 2 + s(5 + 2c 2 , - 3ft 2 - c e {3 + p 2 ))) 
-(-1 + c e )s(5 + 2c 2 - 3/3 t 2 - c e {3 + /3 2 )))) + ^6 2 (c 2 (2m 2 (-ll - 2c + 
-(-1 + c e )s(5 + 3(3 2 + c e (7 + P 2 ))) - d 2 (4m 2 (7 - 2c e + 2c 2 + /3 2 ) 
+(-1 + c g )s{5 + 3(3 2 + c e {7 + (3 2 )))) + ^a 2 Ub 2 m 2 (l + 6c g + 3c 2 - 2(3 2 ) 

-d 2 (2m 2 (3 + 18c e + 3c 2 g - 8(3 2 ) + (-1 + c e )s(5 + 3[3 2 + c e (7 + /3 t 2 ))) 
+c 2 (5s + Am 2 t (3 2 + 3s/? 2 - 2c 9 (24m 2 + s(-l + $)) + c 2 g (12m 2 - s(7 + ft 2 ))) 

1 / TO 2 

2 

(a 3 c + ab 2 c + a 2 c 2 - b 2 d 2 )s + + c e ) 2 (a 2 c 2 - b 2 d 2 )s 2 (5 + 2c e + f3?) 



im 2 (-5 + 8c e + c 2 g -Al3 2 ) 



(Tft \ r 
-^J s ]32(a 4 + b 4 + 2a 3 c + 2ab 2 c + a 2 {2b 2 + c 2 ) - b 2 d 2 )m 2 + 8(-3 + 2c g + c 2 e ) 



h = tt) .s 2 (a 2 + 6 2 ) 2 (-l + Ce ) 2 



(A3) 
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